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Abstract 

This work investigates the consequences of imposing a volume constraint on the maximum power 
'"^ ■ that can be absorbed from progressive regular incident waves by an attenuating line absorber heav- 

ing in a travelling wave mode. Under assumptions of linear theory an equation for the maximum 
absorbed power is derived in terms of two dimensionlcss independent variables representing the 
length and the half-swept volume of the line absorber. The equation gives the well-known result 
for a point absorber wave energy converter in the limit of zero length and it gives Budal's upper 
bound in the limit of zero volume. The equation shows that the maximum power absorbed by a 
^.^i heaving point absorber is limited regardless of its volume, while for a heaving line absorber whose 

(~| I length tends to infinity the maximum power is proportional to its swept volume, with no limit. 

Power limits arise for line absorbers of practical lengths and volumes but they can be multiples of 
those achieved for point absorbers of similar volumes. This conclusion has profound implications 
for the scaling and economics of wave energy converters. 

Keywords: Wave power absorption, attenuating line absorbers, point absorbers, volume 
constraints. 
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1. Introduction 



2. Background Theory 



A large variety of wave energy converters 
(WECs) have been proposed, with a few demon- 
strated at full scale. This paper is concerned 
with the following two categories of buoyancy- 
driven heaving WECs: the archetypal point ab- 
sorber which is a single scmi-submcrgcd float 
with horizontal dimensions a that arc small com- 
pared to the wavelength A of incident wave; and 
the archetypal attenuating line absorber which 
is a continuous line of scmi-submcrgcd floats 
with its long dimension orientated in the dir- 
ection of travel of the incident wave, its width a 
small compared to the wavelength and its length 
I of the order of a wavelength. 

It is well established from linear water-wave 
theory that the upper limit of the capture width 
of a heaving point absorber is \/2tt (see, for 
example, Evans 0, Mei [||, Newman [||). As 
A increases, however, reaching a capture width 
of A/27r requires an increasingly large swept 
volume. Of more practical significance arc the 
results of Evans and Pizer that account 
for a constraint on the motion — and, therefore, 
the swept volume — of the float. 

Farley |^ provides a theory for the capture 
width of an attenuating line absorber oscillat- 
ing in a travelling wave mode, but he does not 
incorporate a motion constraint. Newman [Q 
does consider a motion constraint but his the- 
ory is restricted to modes of motion that have 
uniform temporal phase and so he does not solve 
explicitly for the more practical travelling wave 
mode. The present work follows the approaches 
of Farley |^ and Rainey |^ in considering the 
travelling wave mode response, and it derives 
the maximum capture width with a motion con- 
straint to limit the maximum swept volume of 
the device. 

The method calculates the capture width by 
considering the radiated and difi^racted waves 
from a line composed of a continuum of point 
wave sources and their interaction with the in- 
cident wave in the far-field. The motion con- 
straint is applied by matching the far-field ra- 
diated wave amplitude with the swept volume 
of a line segment using the relative motion hy- 
pothesis. Example calculations are provided for 
point and line absorbers of different volumes in 
wave conditions of practical interest. The effects 
of variations in wavelength and wave height are 
also illustrated. 



2.1. Coordinate systems 

The coordinate systems and related notation 
arc described here. Let S denote an iner- 
tial reference frame described by the Cartesian 
coordinates (x,y,z) = (xi,a;2,X3) with basis 
{ex,ey,ez) = (61,62,6:3) and by the cylindrical 
polar coordinates (r, 9, z) with basis (e^, eg, e^). 
The origins of these coordinate systems coin- 
cide, the z-axis is directed vertically upward, the 
plane of the undisturbed free surface is at z = 
and the fluid bottom is at z = — /i. When 6 = 
the unit vectors and e^; are equal. For any 
vector n, the following notations are equivalent: 
n • V = TT-Q d/dxa = d/dn = 9„, where V is the 
gradient operator and a can take values from 
{1,2,3} or {x,y,z}. 

2.2. Governing equations 

Consider a rigid body floating on the surface 
of an inviscid, irrotational, incompressible fluid 
and interacting with a plane incident wave. All 
motions are assumed to be time-harmonic with 
angular frequency lo. It is assumed that the 
body-induced perturbations are small so that 
the linearised theory of the interaction of water 
waves and structures can be applied and higher- 
order efi^ects can be neglected. 

Under these assumptions the fluid is described 
by cither its real velocity potential, or by 
the time-independent complex amplitude of the 
velocity potential, 0, where 

$ = 5R{^e-*'^*} . (1) 

The velocity potential, (j), must satisfy the 
Laplace equation, 

V^^^O, -h<z<0. (2) 

In the linearised approximation for small amp- 
litude waves and constant atmospheric pressure 
at 2; = 0, the boundary conditions at the sea 
bed and the free surface are 

^ = 0, z = ~h, (3) 

^Lo^ + g^ = 0, z = 0, (4) 
oz 

where g is the acceleration due to gravity. 

The boundary condition on the wetted sur- 
face, denoted by Sb, of a heaving body is ob- 
tained by equating the normal velocity of the 
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body to the normal velocity of the fluid on Sb- 
The result is 



-iuj UzAy on Si 



(5) 



where n is the unit vector in the direction of the 
normal pointing into the body at any point on 
Sb, nz is the component of n in the direction 
of e^, and is the complex amplitude of the 
oscillatory heave motion of the body about its 
mean position. 

2.3. Extraction of wave energy 

The time averaged power, denoted by P, 
transferred from the fluid to the body is ob- 
tained by integrating the time average of the 
fluid pressure, p, times the fluid velocity, d^/dn, 
over Sb. This may be expressed a^ 



P 




(6) 



where T = 2tt/uj is the period of oscillation. 
Note that P is the time average of the rate of 
change of work done on the body so that a pos- 
itive value of P represents energy absorbed by 
the body from the fluid. 

The pressure field in the fluid is given by the 
linearised Bernoulli equation which can be writ- 
ten as 

— +gz+-^—, 7) 
at P P 

where po is the pressure at 2 = 0, g is the ac- 
celeration due to gravity and p is the fluid dens- 
ity. In the absence of a wave it is assumed that 
p = Po is constant at z = 0. 

Substituting $ from and p from into 
the integrand in (^), it can be shown that 



1 



P^dt 
on 



iujp / d4> - d(j) 



dr 



dn 



(8) 



where (j) is the complex conjugate of 4> and use 
has been made of the fact that the integrals of 
the oscillatory terms are zero. Substituting (|[) 
into (O) gives 



P 



ILOp 



Sb 



90 

dn 



dS. (9) 



^Within the linear theory applied here, it is accept- 
able to integrate over the time average of the position of 
the surface of the body as integrating over the instant- 
aneous position of the surface of the body adds at most 
a second-order term to the time average of the absorbed 
power. 



Application of Green's second identity to the 
integral in (H), and using the homogeneous 
boundary conditions (H) and (Q), gives the in- 
tegral over the body's surface in terms of an in- 
tegral over a far-flcld control surface. The result 



Sb 



dn 



Sc 



dn 



dS 



where Sc is a cylindrical control surface whose 
axis is in the direction and which encircles 
the body. Substituting (|l^) into converts the 
integral over the body's surface into an integral 
over the far-fleld control surface, and it follows 
that 



P = - 



ILOp 



Forthcoming algebra is simplifled by substitut- 
ing 

del) idcj) _ ( d(t> 
dn dn dn 



into (11) to giveQ 



P ~ 

2 



Sc 



dn 



dS 



(12) 



This equation expresses the power, P, trans- 
ferred from the fluid to the body in terms of 
the total power entering the control volume 
bounded by the control surface, Sc. 

2.4. Linear decomposition of velocity potential 

According to the linearised theory of the in- 
teraction of water waves with heaving bodies, 
the total velocity potential of the fluid may be 
approximated by 



(j) = Aqipo + Aoipd + Az (fz, 



(13) 



where: (po and (pd are the velocity potentials 
per unit complex amplitude for the incident and 
diffracted waves respectively; Aq is the complex 
amplitude of the incident wave; ipz is the ve- 
locity potential per unit complex amplitude of 
the oscillating heave motion of the rigid-body; 



^Note that (|ll|) and (^) appear asvmnumbered equa- 
tions at the top of page 320 in Mci |H|. 
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and Az is the complex amplitude of the oscillat- 
ing heave-mode of the body as introduced pre- 
viously in (^). Each of the velocity potentials 
(po , ipd and ifz satisfy the Laplace equation 
the sea bed condition (H), and the free-surface 
condition (Q). Furthermore, (pd and ipz satisfy 
the Sommcrfield radiation condition at large dis- 
tances from the body. 

In this treatment the incident wave is a plane 
wave travelling in the positive a;-direction with 
a velocity potential given by 

, A ■ A 9 cosh(fc(z + h)) jt,^ 

a; cosli(Kn,) 

(14) 

and a dispersion relationship given by 

= gk tanh(/c/i), (15) 

where lo is the wave's angular frequency and k 
is its wavenumber. 

The linearised boundary conditions on Sb are 
obtained by substituting from ( |l3| ) into the 
boundary condition (||). Setting = Q gives 



n • Vi^d = -n • V(^o on 5*5, 



whilst setting — gives 
n • V(^z = —iu Tij 



on Sb 



(16) 



(17) 



2.5. Relative motion hypothesis 

This section summarises the application of the 
relative motion hypothesis and matched asymp- 
totic expansions used by Mclver to give the 
leading order relationship between the near and 
far-field expressions of ipd and Lp^. Note that 
the finite-depth analysis of Melvei §9 and 
§10] is not applicable to the problem here as it 
applies to cases in which the body dimensions 
are of the same order as the water depth. We 
therefore proceed with the finite-depth gener- 
alisations of the Mclver 's infinite depth results 
§7]. 

It is assumed that Sb is centred at the ori- 
gin of the coordinate system and that its extent 
is small compared with the wavelength of the 
incident wave. Under these assumptions it is 
consistent with the analysis presented here to 
approximate (po on iSs by using just the terms 
up to first-order in ka in the Taylor expansion of 
(po from (14) with respect to the dimensionless 
variables kx and kz about the origin x = y = Q. 
That is, 



'i— (l + k (ix 

CO 



tanh(A:/i)) 



0{kaf) . 
(18) 



This is the finite-depth equivalent to (158) in 
(l^. Equation (|l|) can be substituted into the 
body boundary conditions (^6|), (p7|), and the 
equivalent surge boundary condition for (/j^, to 
show that in the vicinity of the body (pd may be 
expressed as a linear combination of ipx and ipz , 
that is, 



•fx 



tanh(/s/i) 



LPz on Sb- 



This is the finite depth equivalent to the near- 



field relationship given by (160) in |10 



The far-field relationship between the dif- 
fraction and radiation potentials is obtained 
by Mclver flQ, §7] using matched asymptotic 



analysis. He shows that for an axisymmetrie 
surface-piercing body the far-field velocity po- 
tential of radiated waves arising from heave mo- 
tions is proportional to (fca)^, whereas that from 
surge motions is proportional to (kaf' . In par- 
ticular, 



using 



Melvei's ]10[ far-field expressions 



(132) and (161), it is clear that the leading or- 
der terms in the far-field diffracted and radiated 
waves are related by ipd ~ —fz- Thus, when 
performing the integral over the far-field control 
surface Sc in ( p^ ) it is consistent with the ap- 
proximations used here to write 



on Sr_ 



(19) 



3. Heaving attenuating line absorber 

3.1. Capture width for unlimited volume 

The slender body approximation is used to 
derive an expression for the capture width of 
an attenuating line absorber, as in Newman Q]. 
This approximation gives an expression for the 
far-field radiated velocity potential at a distance 
r from a slender body. For a body of length 
I and beam b it states that, provided 6 <C A 
and / -C r, an approximation for the far-field 
radiated velocity potential at r can be obtained 
by integrating over the contributions from the 
infinitesimal elements of waterplane area of the 
body, b6l where SI ^ A. Denoting by S(t)n,z the 
velocity potential from the heave motions of the 
nth infinitesimal element, it follows that the far- 
field velocity potential due to the heave motions 
of all the infinitesimal elements comprising the 
slender body is given by 



(20) 
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To express each Sipn^z, consider the only 
source of waves in a fluid to be Az amplitude 
heave oscillations of the nth infinitesimal ele- 
ment of the body. Let {r'^,9'^,z') be the polar 
coordinates in the coordinate system for which 
the infinitesimal element is centred at =0. 
The velocity potential per unit complex amp- 
litude of this infinitesimal element at a large dis- 
tance from the centre of the body is given byQ 



1/2 



cosh(fc(z' + /i)) 



cosh(fcft,) 



Q^kr'„-^^/4^ fc^^ ^ (21) 



where 5Hz{0'^) is the Kochin function [for a 
definition, see ^ , which is a dimcnsionless func- 
tion describing the angular dependency of the 
radiated wave amplitude of the heave motion in 



the limit kr 



oo. Generally, for a heaving 



body with infinitesimal waterplane area SSw, 
the Kochin function is independent of 9 and is 
given by SHz = SS^ [see on page 22 of 
[Tot . Thus, for the infinitesimal rigid elements 
comprising a slender body of constant beam, 
SSw = bSl, it follows that 



SHz = -k^bSl. 
2 



(22) 



Here, the assumption is made that the water- 
plane area of the element is not a function of 
the depth of submergence of that element. This 
assumption is strictly valid either for infinites- 
imal relative motions or for vertical wall-sided 
elements or both. 

To apply the slender body approximation in 
( pO| ) , let the coordinates in E of the centre of the 
7ith infinitesimal rigid-body element be x = x„ 
along the line z = y = such that —1/2 < 
Xn < V2 and let the length of the element be 
SI = Sxn- Denote by 

S4>n,z{x, y, z; Xn) = An,z S(pn,zix, y, z; Xn) (23) 

the far-field velocity potential at {x,y,z) in E 
due to heave oscillations of the infinitesimal 
source at a; = a;„ . Let SJ^ denote the coordinate 
system with the same Cartesian basis vectors 
as S but with its origin shifted to (a;,i,0,0) as 



^See, for example, Mclver ^xj], in particular see 
Eqs. (10) , (143) and the unnumbered equation immedi- 
ately following (14.S) w hirh is the Kochin function, de- 
noted by ^3 in Mclver 's notation. 



measured from S. Thus, in E'^ the ?ith infin- 
itesimal element is at (xjj, y^j, z^) = (0,0,0) in 
Cartesian coordinates or, equivalently, at r'^ ~ 
in polar coordinates z'^). The Cartesian 

coordinate transformations between E' and E 



are 



y, 

z. 



(24) 
(25) 
(26) 



Substituting (g4D-(|2^) into r^- 
cxpanding, and using = : 
r cos 9, it can be shown that 



Xn COS 9 + 



y- and a; 



(27) 



Substituting (|22|), ( p7| ) and 61 = 6xn into 
( PH ) and neglecting terms of order 0{x^/r) and 
higher, (5(^„^z(r^, ^^j, z) is transformed from the 
E^ reference frame to the E reference frame 
yielding 



Sipn,zir,9,Z]Xn) ~ -■ 



gk^bSxn 



uj(2T:krfl'^ 

cosh(fc(z + fe)) ^,fc^_,^/4^-.fc:,„ cose 
cosh(A:/i) 

kr 00. (28) 

For an attenuating line absorber composed of 
heaving elements oscillating with relative phases 
that give the dynamics of a travelling wave mov- 
ing in the positive x-direction along the length 
of the attenuator, the complex amplitude can be 
written as the function of .t„ given by 



An,z{Xn) — Az e 



(29) 



where Az is the complex amplitude of the heave- 
mode oscillation of the element located at the 
origin of the E reference frame. Substituting 
(H) into (|23|), it follows that 



S(f>n,z{r, 9, z; x„) = Az e''""'" Sipn,zir, 9, z; a;„). 

(30) 

From (|o|), (|2|) and (l)z{r, 9, z) can be eval- 
uated as 



c^z{r,9,z) 



kl 



gk'^Azbl cosh{k{z + h)) 
uj{2TTkrf'^ cosh(fc/i) 



jo( y (1 - cosO) ] e*'='-W4^ ^ 00, (31) 
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Figure 1: Schematic of the surface elevation of the wave 
radiated from a line of heaving elements oscillating on 
otherwise still water. The line has a length of i/A = 2 
and is oscillating in atravelling wave mode with velocity 
potential given by (pl|). 



where jo is the zeroth-order spherical Bessel 
function of the first kind (also known as the 
sine cardinal or sine function). All of the 
^-dependency in ( ^l| ) is within the spherical 
Bessel function. The heavc-modc oscillations 
have wavelength A and are travelling in the e^:- 
direction along the length, I, of the line ab- 
sorber. Figure |l| shows a schematic of the 9- 
dependency of the surface elevation of the radi- 
ated wave, given by 77 = —g~^ f^0z/'9^L=o' ^r 
(/)^ in (31) with l/X = 2. The concentration of 
waves in the e^^-direction resulting from these 
phased motions is clearly visible in the figure. 

For a given amplitude Az and waterplanc area 
bl, and since jo{0) = 1, it can be seen from ( ^T| ) 
that in the limit kr 00 the ratio of 02 for 
a line absorber of length / to that of a point 
absorber is simply joi^ (1 ~ coa6)). Since the 
surface elevation of the radiated waves is pro- 
portional to 02 at z = 0, the ratio of the surface 
elevation of a line absorber to that of a point 



absorber is also jo{^Y ~ ^'^^^)) ■ Graphs of 
jo(^ (1 — cos 6*)) for a point absorber and four 
line absorbers of different lengths are shown in 
Figure ||. From this figure it is clear that the 
point absorber radiates waves in all directions 
with the same amplitude, whereas the line ab- 
sorbers radiate at the same amplitude as the 
point absorber in the = or e^j-direction, but 
they radiate at reduced amplitudes in other dir- 
ections. Also evident is that the longer the line 
absorber, the narrower the "beam" of radiated 
waves and the less the energy radiated. 

An expression for the full velocity potential, 
(j), resulting from heave motions and wave dif- 
fraction in an incident wave is obtained by sub- 
stituting the relative motion hypothesis ( p^ 
into (O) to give 



• -■ //x=i 



— ■ /A = 4 ' 



TTTT — ^ 

point / ■/■ 

111=0.5 ' I 
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Figure 2: Plots of the function jo(^ (1 - cos9)j. This 
function describes the 6-dependency of the surface eleva- 
tion of the radialfid waves and that of the velocity poten- 
tial, 02, from (Bll). The 9-dependency plots for a point 
absorber and four line absorbers of different lengths I / A 
are shown. 



where the relative complex amplitude is defined 
by 

Ar = Az- Aq. (33) 

Substituting ( |3^ ) into the integrand in the ex- 
pression for absorbed power in (|l^) and expand- 
ing gives 



<i>^ = \Ao\ (^0-^ VA^Ar^Pz-^ 

on an on 

+ AoAr^O^ + \Ar\\z^. (34) 

The incident wave will contribute no net energy 
fiux through the control surface Sc, therefore 
IIsc 'Pod^o/dndS ~ 0. With this observation 
the term proportional to \Aq\^ in ( ^ ) can be 
ignored and the integral in (g^ can be written 



Sc 



90 
dn 



dS = 



Sc 



AoArifO 



difiz 
dn 



4> ^ Ao(po+ Ar ipz 



(32) 



+AoArVz ^ + I A, I' ^z ^) dS. (35) 
on On J 

The terms on the right-hand side of ( |35| ) are 
evaluate d in cylindrical polar coordinates in A p- 
pendix A. The result shows that P in (|12^can 
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be written as 



where Cg is the group velocity given by 



P = |3 { {AoAr - AoAr) Cggk hi 



^i\Ar\^ Cggk'^h^P I{kl)Y (36) 
where I{kl) is given byQ 

I{kl) ^ i- j1 (y (1 - cos 9)^ de, (37) 
4 

= - cos(fc/) Jo(kl) 
3 

2 

+ (2klsm(kl) - cos( kl)) Ji(kl), 

3kl 

(38) 

and where Jq and Ji are the zeroth and first- 
order Bessel functions of the first kind. Since 

AoAr - AoAr = 2i3 {^oA} , 
it follows that 

P = ^Cgpgkbl (2Q{AoAr} - \Ar\^ k^bll{kl)^ . 

. . . 

In this equation the first term, which can be 

positive or negative depending on the relative 
phase of Aq and Ar , gives the energy transferred 
from the fiuid to the body. From the point of 
view of a near-field perspective, it represents the 
fluid pressure on the body times the velocity of 
the body; from the point of view of a far-field 
perspective, it represents the interference of the 
incident and radiated waves in the far-field. The 
second term, which is always negative, gives the 
power lost due to energy radiated away from the 
body. 

The mean power absorbed by the body, P in 
( p9| ) , is equal to the mean power extracted from 
the incident wave. The mean power extracted 
from a wave can be written as the product of 
the mean energy fiux in the wave, denoted by 
J, and the power capture width, denoted by w, 
that is 

P = wJ. (40) 

The mean energy flux of the plane incident wave 
described by (P_4h can be written as 



J ^ T;CgPg\Ao\ 



(41) 



^The analyticaLcxDr 
Tp is given by Farley |198 



in (Bq) for the integral in 
Eq. (20)]. 



duj I I 

^ ~dk^ 2k 



2kh 



sinh(2fc/i)) 



(42) 



Substituting (H) and (|l]) into and rearran- 
ging gives the power capture width as 



2k bl 



w 



\Ao 



\Ar\ 



{AoAr}-kH^l^^Iikl). (43) 

1^0 



This capture width can now be optimised with 
respect to the phase and amplitude of Ar . Since 
the second term on the right-hand side of (^3|) is 
positive and independent of phase, to optimise w 
with respect to phase requires the maximisation 
of 

^{AoAr} ^lAoWArlsmA (44) 

with respect to ipr = arg {AoAr) — arg (^o) ~ 
arg (Ar), which is the relative phase between 
and Ar- Maximisation requires that sinipr = 1, 
implying that ipr = 7r/2 so that Ar lags Ao by 
7r/2 radians and the response velocity (which 
is proportional to dnfpz) is in phase with the 
excitation force (which is proportional to 0o)- 
Substituting t/j^ — -k 12 into (p4[), and the result 



into (43), gives 



I{kl). (45) 



The flrst term in this equation, which is always 
positive, accounts for the net power entering the 
control volume bounded by Sc due to the inter- 
action of the incident and radiated waves; the 
second term, which is always negative, accounts 
for the energy leaving the control volume due to 
the radiated waves. 

To optimise (|5|) with respect to the relative 
amplitude j^^l consider 



d\Ar 



2fc3 6^/2 



\A^ 



the solution of which is 



\Ar\ = 



1 



fc2 hll{kl) 



■l{kl) = 0, 



(46) 



Substituting \Ar\ from ( ^ ) into the first and 
second terms on the right-hand side of (^5|) gives 
the optimal amplitude condition that the radi- 
ated power represented by the second term is 
exactly half of the absorbed power represented 
by the first term. This, together with the op- 
timal phase condition that Ar lags Ao by 7r/2 
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radians corresponds to the familiar condition of 
"impedance matching" to maximise the energy 
transfer between a source and a load. 

The value of \Ar\ in ( p6| ) is not limited by the 
volume of the device, and since k is typically re- 
latively small the value of jA^j can be very large 
in comparison with |^ol- Because of this, when 
\Ar\ from (|4^ ) is substituted back into ( ^5|) it 
gives the optimum capture width for a device 
of unconstrained motion, or equivalently, unlim- 
ited volume, as 



1 



kl{kl) 



if lArl is unlimited. 



(47) 



This is identical to the result reported by Farley 
[H (see Eq. (18) and the appendix of that work). 

In the limit that the length of the line ab- 
sorber tends to zero, by ( A. 12 ) in Appendix A , 
limfe;_s.o I{kl) — 1 and therefore w = 1/k, which 
agrees with the result obtained for the capture 
width of an unconstrained heaving point ab- 
sorbed by, for example, Evans [Q and Newman 

i- 

At this point a physical interpretation can 
be given as to why higher capture widths are 
achievable from line absorbers than from point 
absorbers: for a line absorber the energy being 
carried away by radiated waves can be made ar- 
bitrarily small while maintaining the amplitude, 
\Ar\, of the radiated wave in the 9 = direction 
that is necessary to remove energy from the in- 
cident wave by destructive interference. This in- 
terpretation is clear from consideration of (^Sj) 
in the light of (37) and Figure |[ Obviously, 
to increase the capture width the energy loss 
given by the second term on the right-hand side 
of ( ^5| ) should be minimised while maintaining 
a high value for the energy absorption given by 
the first term on the right-hand side. For a point 
absorber I{kl) = 1, so \Ar\ is the only para- 
meter available that is independent of the incid- 
ent wave and whose value can be adjusted to 
maximise w. The magnitude of \Ar \ needs to be 
increased to increase the energy absorption from 
the first term in (|45|), which is linear in \Ar\] but 
if \Ar\ is made too large the energy loss from 
the second term in (p5|), which is quadratic in 
\Ar\^ dominates. For a line absorber, however, 
the added degree-of-freedom of length can be 
used to reduce the spread of the radiated waves, 
as shown in Figure |^. This reduces the mag- 
nitude of I{kl) while keeping a large magnitude 
for \Ar\, so achieving a high energy absorption 



from the first term on the right-hand side of ( p5| ) 
and a low energy loss from the second. 

3.2. Capture width for limited volume 

In the preceding section no constraints were 
placed on the maximum value of the relat- 
ive amplitude, \Ar\^ and the maximum capture 
width for a line absorber of unlimited volume 
was derived as (^). In this section a motion 
constraint is placed on j^^l in (|45|). In the lin- 
ear theory presented here such a constraint on 
\Ar \ may be equated to a limit on the maximum 
swept volume of the device. This constraint 
may be a consequence of the limited volume of 
a device, or it may be a consequence of other 
engineering or control constraints that limit the 
range of heave motion. Assume for simplicity 
that the depth, D, and the beam, b, of each 
vertically-wall-sided element are independent of 
its position along the length, I, of the line ab- 
sorber. Here, D is the vertical difference in 
position between the minimum and maximum 
submergences of each element of the line ab- 
sorber. If the element is symmetric about the 
plane z = the draft and the freeboard are equal 
and their sum is equal to D. In this case half 
the maximum swept volume of the line absorber 
is given by 

v^4, 



and the maximum value of l^^l is 



max (|Ar|) 



D _V 

'2 ^Vr 



(48) 



Substituting \Ar\ from (^ into (]4^) gives the 
volume limiting condition in terms of \Aq\ as 



■.{\A„\)^k^VI{kl). 



(49) 



Thus, for a limited volume device the maximum 
\Ar\ is reached when \Ao\ > k^V I{kl) and the 
maximum capture width is given by substituting 
the maximum of \Ar\ from ( |48| ) into (^). This 
gives 

kV / k^V , ■ 

if \Ao\ > k^VI{kl). 

When \Ao\ < k'^V I{kl) the optimum capture 
width is still given by (^7|). Thus, the maximum 
capture width of a volume-limited heaving line 
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attenuator can be written as 
1 



w = < 



ki{kiy 



\Aq\ < k^VI{kl), 



kV / k'^V \ 
where an analytic expression for likl) is given 

by m. 



4. Dimensionless capture width and ab- 
sorbed power 

It is instructive to formulate the results of §|| 
in dimensionless units. The dimensionless inde- 
pendent variables of device length and device 
half-swept volume, and the dimensionless de- 
pendent variable of capture width, are defined 

by 



r 

V* 



ki, 

k^V 
kw. 



(51) 
(52) 
(53) 



Since capture width and absorbed power are re- 
lated by (^, the definition of w* in (^ ) sug- 
gests a non-dimensional power defined by 



p* 



kP 



(54) 



where J is given by (p]). From (^), (^^ and 
( ^ ) is follows that the dimensionless variables 
w* and P* are equivalent. Substituting (^l])- 
( ^ ) into (^o|) gives the dimensionless capture 
width and power as 



1 



= P* = 



i{i*y 



V* i{r) > 1, 



y*(2-y* /(/*)), T/*/(r)<i, 

(55) 

where / (/*) — I{kl) is the integral given by (^). 
The forms of the functions I {1*) and l/I {1*) are 
shown in Figure H. In these plots I* j^-n ~ l/X is 
used as the dimensionless unit on the abscissa 
as these units correspond to multiples of the 
wavelength. Shown in the right-hand side plot 
in Figure H is the line V* = l/I (l*) which marks 
the boundary between the limited and unlimited 
volume regimes. 



The point absorber limit is obtained by sub- 
stituting limi._j.o ^ (^*) ^ 1 into (p5|) to give 



lim w* 

1*^0 



1, 



2V*~V* 



V* > 1, 
V* < 1. 



(56) 



The term in ( p6[ ) which applies when V* >1, 
that is w* = 1, is the well known result for the 
theoretical maximum capture width of a heav- 
ing point absorber of unlimited volume. When 
V* 1 the term, which is first-order in volume, 
dominates and the capture width can be writ- 
ten as w* = 2V* , or, equivalently, P* /V* — 2. 
This corresponds to Budal's upper bound^ on 
the maximum capture width, or power per unit 
volume, of a finite-volume heaving point ab- 
sorber [see |ll|. The term that is second-order 
in volume is a modification arising from the in- 
clusion of energy loss due to radiated waves. 

The infinite length line absorber limit is ob- 
tained by substituting lim;.^oo ^ ('*) — > into 
(B3) to give 



lim w* 2V* 

!* — >oo 



for ah V*. 



(57) 



This can also be written as P* /V* = 2 for all 
V*. Comparing ( p6| ) and (|^ shows that for 
a heaving point absorber the maximum cap- 
ture width is limited to ui* — 1, no matter 
how great the volume; for a line absorber in 
the limit infinite length, however, the maximum 
capture width scales linearly with volume, that 
is w* ^ 2V*. 

Plots of dimensionless capture width (or ab- 
sorbed power) from (^5|) are shown in Figure ^. 
The two plots give different illustrations of the 
effect on w* of changing the volume and length 
of a device. The solid line in the left-hand plot 
in Figure ^ shows that increasing volume above 
V* = 1 gives no increase in w* for a point 
absorber; the dashed lines show that, for line 
absorbers, increasing V* above 1 can increase 
w* , and that the increase is greater for a longer 
device. In fact, the dotted line marks where 
w* = 2V* which shows that the maximum di- 
mensionless capture widths for an infinite length 
line absorber is equal to its dimensionless swept 



^Budal's result is generally quoted as the upper 
bound on the absorbed power per unit volume given by 



-pgui \Ao \ , 



P 1 

where Vt = 21/ is the total swept volume of the body. 
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Figure 3: The dimensionless functions I (kl) = I (l*) and V* = 1/ 1 {I*)- 



volume. The fine dotted line marks where w* = 
V* which is the line of transition between the 
unlimited volume regime below the line and the 
limited volume regime above it. It is clear from 
this plot that the minimum half-swept volume 
required for a point absorber to achieve its max- 
imum capture width is V* = 1, and for a line ab- 
sorber it is V* = w*. Note that V^* = 1 is equi- 
valent to V/X^ \ Ao\ = l/47r2 w 0.025, which aids 
in interpreting the dimensionless volume as it 
states that the minimum half-swept volume for 
a point absorber to achieve its maximum cap- 
ture width is about 2.5% of the natural measure 
of volume given by the square of the wavelength 
times the amplitude of the wave, that is, |^o|- 

The right-hand plot in Figure ^ shows the 
gains in w* that can be achieved by "stretch- 
ing out" different fixed-volume point absorbers 
to line absorbers of different lengths. The points 
where the lines terminate at Z* =0 represent the 
dimensionless capture widths of point absorbers 
of the appropriate volumes. There is little dif- 
ference in w* between point and line absorbers 
when V* < 0.5, that is, well within the volume- 
limited regime of the point absorber. For the 
unlimited volume regime of 1^* > 1 the capture 
width of the point absorber is at its theoretical 
maximum of = 1 and the lines for V* = 1, 
2 and 4 all converge to w* = 1 as — > 0. For 
finite length line absorbers, however, significant 
increases in capture width are achievable above 
w* = 1, as shown in the figure and in the ex- 



half-swept 


length. 


capture- 


volume, V* 


r/27r = l/X 


width, w* 


1 





1 


1 


1 


1.684 


1 


2 


1.782 


2 





1 


2 


1 


2.735 


2 


2 


3.127 


3 





1 


3 


1 


3.154 


3 


2 


4.036 


unlimited 





1 


unlimited 


1 


3.162 


unlimited 


2 


4.583 



Tabic 1: Examples of maximum dimensionless capture 
widths for different dimensionless volumes and lengths. 



amples given in Table § 

Plots of dimensionless absorbed powers (or 
capture widths) per unit dimensionless volume 
are shown in Figure |5[ In the left-hand plot 
in Figure ^ the fine dotted line marks the line 
of P*/V* = 1, which is the line of trans- 
ition between the limited volume regime above 
the line and unlimited volume regime below 
it. In this plot the limited- volume dimension- 
less powers converge to P* /V* = 2 a.t V* = 
for all device lengths I*. This maximum 
value of P* /V* = 2 corresponds to Budal's up- 
per bound expressed in the dimensional units 
of (|l|)-(|l. Budal's upper bound applies to 
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V = 025 
V' = 0.5 

v = -\ 

V = Z 

V = 4 

Unlimited voi. 



1 

2 

/•/2it 



Figure 4: Dimensionlcss capture width, 
dimensionless device length, l* . 



w* , from fcSh as a function of dimensionless half-swept volume, V, and 



all wave energy absorbers, however, point ab- 
sorbers can only reach Budal's upper bound in 
the zero- volume limit, whereas line absorber in 
the infinite-length limit can reach Budal's upper 
bound for any volume. This is evident from the 
infinite- length limit of ( p5| ) giving P* /V* = 2 
which is independent of volume. From the left- 
hand plot in Figure || it is clear that for small 
volume devices there is no significant advantage 
to using an attenuating line absorber compared 
to a point absorber. For larger volume devices, 
however, the same plot shows how the power 
per unit volume decreases less rapidly for a line 
absorber than for an equivalent volume point 
absorber, and the rate of decrease of power per 
unit volume is less for longer line absorbers (be- 
ing zero in the limit of an infinitely long line 
absorber). In fact, by dividing volume of a 
single point absorber with V* = 1 into ever 
smaller volume point absorbers the maximum 
power per unit volume that can be achieved goes 
from P*/V* = 1 to P*/V* = 2 and, therefore, 
the maximum absorbed power is doubled. This 
severe drop-off in P* /V* with increasing V* for 
point absorbers is an illustration of the "small is 
beautiful" argument made by Falnes . 

The right-hand plot in Figure |^ gives an 
alternative illustration of how, if efficiency is 
measured in terms of absorbed power per unit 
volume, for all device lengths, 1* , devices of 
smaller volume always have greater efficiencies. 
Note that in this plot the dotted line starting at 



half-swept 






volume. 


length. 




V (m^) 


I (m) 


WEC of similar size 


240 





OPT PB150 |12 


940 





WaveBob |13| 


350 


120 


Pelamis FSP JlJ 


790 


180 


Pclamis P2 ||l5[ 


1700 


210 


possible future Pelamis 



Table 2: Selection of WEC sizes chosen to represent real- 
istic production sizes based on estimates of the sizes of 
devices already built or recently proposed. 



/* = and P* /V* = 2 represents the power per 
unit volume in the unlimited volume regime for 
the case V* = 0.5. The corresponding line for 
V* = 0.25 is off the scale. 



5. Application 
WECs 



to realistically sized 



The advantage of formulating the results in 
dimensionless units is that the number of in- 
dependent variables is reduced from four (fc, 
V, and \Ao\) to just two (/* and V*). In di- 
mensionless units, however, it is not immedi- 
ately obvious where a particular incident wave 
and WEC combination will occur on the graphs. 
In this section we consider the five realistic- 
ally sized WECs described in Table ||. It is 
stressed that the results in the following fig- 
ures are not representative of the actual capture 
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Figure 5: Dimcnsionlcss power per unit dimensionless 
half-swept volume, V*, and dimensionless device length, I 

widths or absorbed powers of these devices, but 
only that devices of these approximate volumes 
and lengths have already been built or may be 
built in the near future. 

Figure 1^ shows the points of maximum dimen- 
sionless capture width plotted against dimen- 
sionless volume for the sizes of WECs given in 
Table 1^ operating in the four incident waves with 
(iJ,r) equal to (2m, 8s), (2m, 10s), (3m, 8s) and 
(3m, 10s). These four incident waves broadly 
correspond to those that would occur close to 
the centres of the available energy resource when 
weighted by annual average occurrence. They 
have incident power in ratios of 1, 1.25, 2.25 
and 2.81 respectively. Also, shown in Figure |6| 
are contours of constant dimensionless length. 
The set of four plots show how the dimension- 
less volume and maximum capture width of a 
particular WEC depends on the wave height and 
period of the incident wave. The plots show that 
the largest volume point absorber appears below 
the line w* = V* and is, therefore, operating 
in its unlimited volume regime and unnecessar- 
ily large for all but the (3m, 10s) incident wave. 
The largest volume line-absorber, however, is 
only unnecessarily large for the (2m, 8s) incident 
wave, even though it is nearly twice the volume 
of the largest point absorber. The smallest point 
absorber has its maximum capture width lim- 
ited by its small volume in all but the (2m, 8s) 
incident wave. The mid-sized line absorber has 
less volume than the large point absorber, yet 
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volume, P* jV* , from (Bfj) as a function of dimensionless 



its maximum capture width is about four times 
larger in the (2m, 8s) wave, three times larger in 
the (3m, 8s) wave, twice as large in the (2m, 10s) 
wave, and one and half times as large in the 
(3m, 10s) wave. For further clarity. Figure ^ 
shows points of dimensional maximum absorbed 
power plotted against dimensional volume for 
the size of WECs given in Table || operating in 
the four incident waves as used in Figure ^. A 
set of conveniently spaced contours of constant 
length are also shown in these plots. This set 
of four plots shows how the maximum absorbed 
power of a particular size of WEC depends on 
the wave height and period of the incident wave. 

6. Conclusions 

The theoretical maximum absorbed power 
and capture width of a limited volume attenu- 
ating line absorber heaving in a travelling wave 
mode in the presence of progressive regular in- 
cident wave has been derived in the frequency 
domain using the linearised theory of the inter- 
action of water waves and structures. The res- 
ults are presented in dimensionless form which 
has the advantage of reducing the number of de- 
pendent variables from four to just two: dimen- 
sionless length, /*, and dimensionless half-swept 
volume, V* . In the zero-length limit the res- 
ults for the limited volume line absorber reduce 
to those for a point absorber; and in the zero- 
volume limit they reduce to the result expressed 
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Figure 6; Dimensionless capture width, w* , from (^) as a function of dimensionless half-swept volume, V* , 
WECs in Table ^ operating in four different incident wave defined by their T). 



for 



by Budal's upper bound. It is shown that, in 
dimensionless units, the maximum dimension- 
less capture width of a point absorber is = 1 
and that the smallest volume required to achieve 
this capture width is y* = 1. The dimensionless 
capture width and volume both depend on the 
incident wave: a value of w* = 1 corresponds 
to a width of about 16% of the wavelength; a 
value of y* = 1 corresponds to about 2.5% of 
the rectangular volume given by the square of 
the wavelength times the wave amplitude. In- 
creasing the volume of a point absorber beyond 
V* = 1 gives no increase in capture width. For 
an attenuating line absorber, however, in the 



limit of infinite length the maximum capture 
width is w* =■ 2V*. Thus, there is no hmit 
to the capture width of a line absorber provided 
it has sufficient volume and length. It is shown 
that a particular fixed-length line absorber will 
have a maximum capture width of w* > 1 and 
that the smallest volume required to achieve this 
capture width is V* ~ w* . 

This has profound implications for the eco- 
nomics of power generation from wave energy 
converters. Even though small volume point 
absorbers are more efficient than larger volume 
point absorbers in terms of power absorbed per 
unit swept- volume of device, engineering limit- 



13 



H = 2m, 7 = 8s 



H = 2m, T= 10s 



V{m'),l{m) 

240, 
t£> - A 940, 

• 350, 120 

1 790, 180 
m - ■ 1700,210 



-A- 



500 



point 

; = 50m 

•-■ ;=100m 

; = 200 m 

; = 400m 



1000 1500 
V(m^) 



I 

1000 



H = 3m, T = 8s 



H = 3m, T= 10s 




— I — 

1000 




Figure 7: Absorbed power, P, from as a function of lialf-swcpt volume, V, for WECs with sizes described in 
Table ti operating in four different {H, T) incident waves. 



ations make installation and operation of very 
large numbers of very small devices excessively 
costly and uneconomic. A point absorber has a 
limit on its capture width and usable volume, 
whereas, in theory, line absorbers can be in- 
definitely scalcd-up in volume and length to 
give unlimited capture widths. This theoret- 
ical advantage of line absorbers over point ab- 
sorbers applies to realistically sized wave en- 
ergy converters currently in operation. Thus, 
line absorbers can be progressively scaled-up 
in volume and length to give increasingly large 
capture widths while retaining dimensions com- 
patible with cost-effective engineering. For ex- 



ample, a line absorber with a volume equal to 
the maximum usable volume of a point absorber 
and a length equal to twice the wavelength of 
the incident wave can absorb nearly 80% more 
power than a point absorber of the same volume. 
Doubling this line absorber's volume gives over 
three times the absorbed power of the point ab- 
sorber, and tripling it gives over four times the 
absorbed power. Thus, it is clear that by in- 
stalling fewer larger line absorbers valuable eco- 
nomies of scale can be achieved from line ab- 
sorbers that are not achievable from point ab- 
sorbers. 
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Appendix A. Integrals by the method of 
stationary phase 



The terms on the right-hand side of (^5|) are 
evaluated in cyhndrical polar coordinates as fol- 
lows. Substituting x = rcos6' into (p^, (po is 
written as 

. g C0Sh(fc(z -t- /l)) ,!.rrnsfl /A 1N 

bj coshykh) 

Also, using dl^ ) along with Ux = cos 9 and Uz = 
on the control surface Sc , and the observation 
that (po is not a function of y, it follows that 



d(po gk cosh{k{z + h)) n ^-tkr cose 



3n w cosh(A;/i) 



on ^c. (A.2) 



Next, using 



n = e^ on Sc, 



and the gradient operator in cylindrical coordin- 
ates 

ar r oO oz 

it follows that 



n-V(p^ = ^— = — — on Ae- 
on or 



can now be simplified and written as 



g^k^bl ( cosh(fc(z -|- h)) 



dn a;2 (27rA:r)^/^ V cosh(fc/i) 



Jo( ^ (1 -cos( 



— iA;r(l — COS 0)+?7r/4 



(A.4) 



dipo 
dn 



g^t^bl /cosh(fc(z + /i)) 



w2 (27rfcr)^/^ V cosh(fc/i) 



d(pz 
dn 



C0S9 Jo( ^ (1 - COS0)^ g.fcr(l-cose)-W4^ 

(A.5) 

2 



ik f gk^hl cosh(fc(z + h)) 



2Trkr 



uj cosh(fc/i) 



Jo( ^(l-cos( 



(A.6) 



In cylindrical polar coordinates the integral 
over Sc on the left-hand of (|3^) is expressed as 



Each of ( A.4)-(A.6) has the same z-dependency. 
Integrating out this z-dependency gives, after 
simplification and substitution of the dispersion 
relationship from (Ra), 



° cosh^(fc(z + h)) ^ _ U!C. 
-h cosh^(fc/i) 



dz 



^ 

gk ' 



(Ai 



where Cg is the gro up ve lo city given by (^). The 
integrals over 9 of (A.4)~(A.6) are evaluated us- 
ing the method of stationary phase, which states 
that for real-valued smooth functions / and g, 
when f'{9o) = and K — !> oo. 



Evaluating d(pz/dn by differentiating ( pl| ) with 
respect to r, it can be shown that 

d(pz igk^bl cosh(fc(2; -|- h)) 

dn w(27rA;r)^/^ cosh{kh) 

jo(|(l-cos0)) e-'='^+-/4 on 5c, 

(A3) 

where jo is the zeroth spherical Bessel function 
of the first kind and, since fcr — > oo on Sc, only 
the highest order terms in kr have been retained. 



Combining @, (Q, ( |A2D and ( |A3D , the 
three integrands on the right-hand side of (SHI) 



g{9)e'''f^'^ de 



27r 



, Kf"{9o) 
with f'{9o) = 0, K ^ oo. (A.9) 

Substituting K = kr, g{x) = jo{ j- (1 - cos 6*)), 
/(x) = cos 6 — 1 and 9o = into (|A.9D gives the 



integral of (A.4) over 9 as 



£.;o(|(i-cos. 



g-ifcr(l-cose) 



-i7r/4 



(AlO) 
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Substituting K = kr, g{x) = 
cos 61 (1 - cos 6>)) f{x) = 1 - cosfl, 
and 9o = into ( |A.9[ ) gives the integral over 9 
of (lO) as 



J cosfljo^yCl 



cos 




(A.ll) 



For convenience of notation the integral of (A.6) 



over 9 is denoted by I(kl) as defin ed by the in- 
tegral in @. As stated by |Farley| [|, Eq. (20)], 
I{kl) can be expressed analytically as in (|38| ) in 
terms of first-order Bessel functions of the first 
kind, Jo and Ji. Since the kl limits of the 
Bessel functions are 

lim Jo(kl) = 1, 



Jijkl) 

kl 



the kl 



lim 



limit of 

lim I(kl) 

kl^O 



(A.12) 



Equations (^), ( |A 10|) ( |A.ll| ) and @ are 
now used, along with (A. 7), to evaluate the in- 
tegrals of ( A.4 )-( |A^ ) over 5c- The results are 



-h 



(fo — — r d9 dz ' 
on 

ipz -T, — r a9 dz ' 
on 

ipz — r d9 dz ' 
On 



Cggkbl 
Cggkbl 



iCggk^b'^P 



■I{kl). 



Substituting these expressions into ( |35| ) and the 
result into (^2|) gives (|36|). 
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